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THE CAUCHY-DIRICHLET PROBLEM 
FOR A GENERAL CLASS OF PARABOLIC EQUATIONS 


PAOLO BARONI AND CASIMIR LINDFORS 


Abstract. We prove regulaiity results such as interior Lipschitz regularity and boundaiy 
continuity for the Cauchy-Dirichlet problem associated to a class of parabolic equations 
inspired by the evolutionary p-Laplacian, but extending it at a wide scale. We employ a 
regularization technique of viscosity-type that we find interesting in itself. 


1. Introduction 


The aim of this paper is the study of the behaviour of solutions to a wide class of 
nonlinear parabolic equations modeled after 


Ut — div - ^Du] =0 in fix := x (0,r) C 

V \Du\ / 


( 1 . 1 ) 


n > 2, where is a bounded domain with boundary and g : ]R_|_ —>■ ]R_|_ is a 
function satisfying 


PO - 1 < Og{s) 


sg'js) 

9{s) 


< Si - 1 


for every s > 0 


( 1 . 2 ) 


with I < go < gi < oo. Notice that we can assume go < gi without loss of generality. 
Indeed, if Og{s) is constant, say Og(s) = p—1 for somep > 1, a simple integration shows 
that g(s) = up to a constant factor, and therefore in this case (1.1) gives back the 
evolutionary p-Laplacian widely studied in particular by DiBenedetto, see the monograph 
[14]. This reveals that (1.1) is a natural generalization of the p-Laplacian, and in effect 
this class of growth conditions was mathematically introduced exactly in these terms by 
Lieberman in [28], even if this kind of condition appears earlier in the applications, see the 
forthcoming lines. 

We stress that quite a comprehensive study of non-negative solutions to the equation 


Uf — div [ip'{u)Du\ = 0 
where the function (p : [0, oo) —>■ [0, oo) satisfies 

SLP^(s') 1 

0 < a < Oin(s) := —< -, for s > 0; 1 + a < Ou,(s) for s > sq 

(p(s) a 


(1.3) 


(1.4) 


for some a S (0,1) and some sq > 0 has been provided by Dahlberg and Kenig [10, 11]; 
see also the books [12, 35]. Clearly, while (1.3) is a generalization of the porous medium 
equation that happens when ip{u) = rt™, m > 0, in the same spirit (1.1) can be seen as a 
generalization of the p-Laplacian. 


As in (1.4), we shall also consider a more stringent growth assumption for g for large 
values of its argument. In addition to (1.2), we shall assume that there exist constants 
ce,e > 0 such that 

g[s) > CIS for any s > 1. (1.5) 
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Note that in the p-Laplacian case (1.5) reads precisely as p > 2n/(n + 2), a completely 
natural assumption in the theory of the evolutionary p-Laplacian operator, see [14, 25, 1]. 
Note moreover that (1.5) is implied by assuming po > 2n/(n + 2), see Paragraph 2.2. 

The regularity for the elliptic and variational counterpart of (1.1) is quite well under¬ 
stood, see for instance [28, 3, 9, 16] for the first argument and [17, 20, 7, 8 ] for the second, 
just to cite some indicative references. In the parabolic setting, however, very few re¬ 
sults are available, and some of them only in particular cases: to our knowledge, only 
[22, 23, 29, 30], all by Lieberman and Hwang, and the recent [5]. 

The difficulty, in particular in finding zero-order results, stems from several facts, the 
main one perhaps being that the equation has very different behaviour, already in the p- 
Laplacian case, in the degenerate (p > 2) and singular (p < 2) cases. In the degenerate case 
phenomena such as expansion of positivity occur, see [ 15, 26], and the diffusion dominates 
[13]. On the other hand, in the singular case the evolutionary character dominates [ 6 ] and 
extinction of positive solutions in finite time could happen, see [14]. In our general setting 
the degenerate case occurs when s i— g{s)/s is increasing, and when it is decreasing we 
have the singular case. However, it might also happen that s g{s) /s has no monotonicity 
whatsoever, making the handling of the equation all the more difficult. The comprehension 
of the interaction of these different phenomena is the key for a better understanding of the 
behaviour of local solutions to ( 1 . 1 ), and in this paper we hope to start to clarify this 
difficult point, which will be the object of future investigations. 

The class of differential operators we study, besides being quite a general extension of 
a well-known operator, finds important applications in the applied sciences, also in view 
of the following observation. Take the convex primitive G of g and consider the general 
minimization problem 



it is often convenient to have energies with a precise dependence on \Du\ of more general 
type than monomial (that is, the case of the p-Dirichlet energy or appropriate extensions). 
For instance, in mechanics, fluid dynamics and magnetism, as first approximation it is cus¬ 
tomary to have dependencies of the energy on the modulus of the gradient of monomial 
type but with exponent depending on the size of \Du\, in order to have mathematical mod¬ 
els fitting the experimental data. In this case g is given by the gluing of different monomials 
(see the example in Paragraph 2.3). At this point, elliptic and parabolic equations having 
the growth described in (1.1) arise naturally as Euler equations or flows of the functional 
in (1.6). In [34], for instance, the two-dimensional stationary, irrotational subsonic flow of 
a compressible fluid is described using an energy defined in the following way: 



for small s, G(s) = quadratic otherwise, (1.7) 


where 7 € ( 1 , 2 ) is the exponent in the law p k, p'^ characteristic of poly tropic gases. 

More in general, see [4, 18, 19], one is lead to consider quasilinear static equations in 
dimension two and three of the type 


div '[p{\Du'^)D'u\ = 0, 


with Du representing the velocity field of the fiow and q = \Du\ being the speed of the 
flow. In this context one introduces the Mach number 


^ [M{q)f := 


(note that we must have p' < 0). In our context, where g{s) = p{s^)s, we compute 
Og(s) = 1 — M{sY. The general theory asserts that a point is elliptic if M < 1 and in 
this case the flow is subsonic, while if M > 1 the point is hyperbolic and the flow there is 
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supersonic. If M = 1 the flow is called sonic. A solution of the boundary value problem 
is called a subsonic (supersonic) flow according to whether all points are subsonic (super¬ 
sonic); note that mixed, or transonic flows can exists, with obvious meaning. However, 
if for some reason we know that the flow maintains a controlled, small speed q, then the 
problem falls in the class of operators we consider; the approximation in (1.7) is a way to 
study flows in the subsonic regime. 


The object of our study will be the Cauchy-Dirichlet problem 


' ut — divA{Du) = 0 

u = '4i 


in Ht, 
on dpflx, 


( 1 . 8 ) 


where A : K” —M" is a vector held modeled after the one appearing in (1.1). In 
particular, we assume it satisfies the following ellipticity and growth conditions; 


(DAiOXA) > 


lei 


\DA{0\ < L 


g(l^l) 

lei 


(1.9) 


for any ^ € R” \ {0}, A € K” and with structural constants 0 < u < 1 < L; we assume 
without loss of generality that Al(0) = 0. The function g is a function as in (1.1), 
satisfying only (1.2) and (1.5). For what concerns tj}, we assume it to be continuous in 
dp^T with modulus of continuity oj^p with respect to the natural distance distpar,G. that is, 
there exists a continuous, concave function : R+ —R+ with uj.pj{Q) = 0 such that 

\ip{x,t) - V'Cy, s)| < u;^(max{|x - y\, - s|)]"^}) 

for every {x, t), (y, s) G dp^lr- As already mentioned, H is a bounded domain of R”, 
n > 2 , whose boundary is of class for some G ( 0 , 1 ); we shall provide some more 
details at the beginning of Section 2. 


In this setting, we state the main result of our paper, which concerns at the same time 
the existence and regularity of a (unique) solution to (1.8). 

Theorem 1.1. There exists a unique solution u, in the sense of Deflnition 2, to the Cauchy- 
Dirichlet problem (1.8), where the vector field A satisfies the assumptions (1.9), with g G 
Ci(R+) satisfying (1.2) and (1.5). In particular, u is continuous up to the boundary and 
moreover if the boundary datum ip is Holder continuous with respect to the natural metric 
distpar,G defined in (2.1), then so is u. 

The following theorem gives some properties together with quantitative estimates for 
the solution described in the previous statement. 

Theorem 1.2. Let u be the solution to (1.8) given by Theorem 1.1. Then u is locally 
Lipschitz continuous and the following estimate holds: 

( r \ 

\\Du\\L^(Q^)<c[f [G{\Du\) + l]dxdt] (1.10) 

V7Q2r / 

for every parabolic cylinder Q2R <£ ^t- The constant c depends on n,gQ,gi,v,L,e 
and C(. Moreover, there exists a modulus of continuity '■ K+ R-p depending on 
n,go,gi,i2, L, e,Q, \\ip\\L<^,ujpj,dD, such that 

\u{x,t) - u{y,s) \ < w„(max{|x - y\, [G"^ (l/|f - s|)] ^}) 

for every (x, t), {y, s) G TTp’. 


( 1 . 11 ) 
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We refer the reader to Paragraph 2.1 for the definitions of the standard parabolic cylin¬ 
ders Qii(xo, fo) and of the parabolic closure of fir- We also mention that in the standard 
case of the evolutionary p-Laplacian our estimate (1.10) gives back exactly the gradient 
sup-estimate available for degenerate and singular equations, see [14, Chapter VIII, Theo¬ 
rems 5.1 & 5.2]. 

Remark 1. Theorems 1.1 and 1.2 hold for a wider class of operators generalizing (1.1), 
which allow the presence of a function g that is not but merely Lipschitz. Indeed, 
we may consider Lipschitz functions g : R_|_ —>■ ]R_|. satisfying (1.2) almost everywhere 
and vector fields A : K" —M" in satisfying the monotonicity and Lipschitz 

assumptions 



( 1 . 12 ) 


for every , ^2 S M" such that |^i | -f |^ 2 1 ^ 0 and for some 0 < < 1 < L. For a proof 

of this fact see the end of Section 6. 

1.1. Novelties and technical tools. We believe that the main interest of this paper, apart 
from the results of Theorems 1.1 and 1.2 themselves (that will be used for instance in 
[31]), is the development of some tools for the treatment of the difficult equation (1.8) 
(see Paragraph 2.3). We prove the Lipschitz estimate as an a priori estimate for problems 
enjoying further regularity. Instead of using a regularization of the type used in [28, 22, 29], 
the regularization we employ is of viscosity type, closer to that in [2]: we consider a vector 
field of the type 


:= (4 * Am + e(l + \^\y e e K”, ^ e (0,1), 


where p ^ 1 is a large exponent and {(pe} a family of mollifiers. This allows us to 
overcome the difficulties of deriving regularity estimates for the approximant problems, 
which we were not able to find in the literature. At this point continuity up to the boundary 
becomes an essential ingredient in the proof of the convergence, as well as the fact that we 
are solving a Cauchy-Dirichlet problem and therefore have a uniform bound on Urtellioo 
given by the maximum principle. 

We use the a priori Lipschitz continuity (and the further regularity) of the approximating 
solutions in a way inspired by [27]. First, we employ the fact that the function v = \Du\'^ 
is a subsolution to a similar problem, see Lemma 3.1. Then, we define an appropriate 
intrinsic geometry (see (3.7)) depending on the growth of the approximating vector field 
Ae, which allows us to rebalance estimates, in the sense that the weight appearing in the 
Caccioppoli estimate for the equation satisfied by v turns out to be essentially constant, see 
(3.8). Here the fact that we can bound the supremum of Du, and thus of v, from above 
is essential. Finally, we conclude the proof using an argument based on an alternative in 
order to get rid of the possible dependence on e in terms of the aforementioned geometry, 
depending in turn on the growth of Ae- 

2 . Preliminary material; notation, the function g, miscellanea 

For what concerns dft, we assume that there exists a radius Rq > 0 such that for every 
point xo G OH there is a unit vector e^g such that the restriction of OH is a graph of a 
function in along the Bxg direction, in the following sense; with T being an 
orthogonal transformation that maps e^o into (0, 0 ,..., 0,1), for every 0 < r < Rq it 
holds 


T ^(OH — Xo) n X (—r, r)) = graphs 
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(see below for the precise meaning of these symbols) with 9 = G 9{B'^) C 

(—r, r) and the norm of 9 uniformly bounded; 

[9]c^,fi < 0 . 

Note that without loss of generality, we can take exg as the inner normal vector in xq- 
{v : {v,exg) = 0} is the tangent hyperplane to Vi in xq', therefore D9{0) = 0. D9 is 
the full gradient of 9 with respect to its n — 1 variables. Finally, by saying that a constant 
depends on 917, we shall mean it depends on 0. 


2.1. Notation. We denote by c a general constant always larger than or equal to one, 
possibly varying from line to line; relevant dependencies on parameters will be emphasized 
using parentheses, i.e., Ci = Ci{n,p, q) means that ci depends on n,p, q. For the ease of 
notation, we shall also use the following abbreviation: 

data := {n,go,gi,i^,L}. 


We denote by 

Br{xo) := {a; € R" : \x — Xq\ < R} 

the open ball with center xq and radius i? > 0; when clear from the context or otherwise 
not important, we shall omit denoting the center as follows: Br = Br{xq). The standard 
parabolic cylinder is defined as 

QRixoRo) ■= Br{xq) X {to — R^,to), 


while we define the natural cylinder as 

Qr{xoM) ■= Br{xo) X {to - 


The latter is strictly linked to the scaling of the equation, see Paragraph 2.6. Unless other¬ 
wise explicitly stated, different balls and cylinders in the same context will have the same 
center. We shall denote, for a factor a > 0, by uBr the ball BaR and by aQR{xo, to) the 
cylinder i?Qi{(xo) x {to — {aR)^,to)', similarly for q;(5^(xo, ^o)- The parabolic boundary 
of a cylindrical domain /C = 22 x F, where 22 is an open domain and F an open interval of 
the real line, is defined as 


dpK. := (22 X infF) U (922 x F). 

Naturally, the parabolic closure of /C is then 1^ := /Cu9p/C. Accordingly with the custom¬ 
ary use in the parabolic setting, when considering a sub-cylinder K. (as above) compactly 
contained in VIr, we shall mean that 22 (e (7 and 0 < inf F < sup F < T; we will write in 
this case K. d Q.t- By 9f7 — xo we mean the set {x € R" : x + xo G 90}. The standard 
parabolic distance is 


distpar((x,7), {y,s)) := max{|x - y\, ^J\t - s|} 

for any (x, t), {y, s) G while a distance strictly related to the scaling properties of 

the differential operator is 


tlistpar.G ((^7 ; (^7 '^)) 


:= max 




( 2 . 1 ) 


Note that (5^(xo, 7o) = {(a:, 7) G R”+^ : distpar,G((a;, 7), (a:o, 7o)) < R,t < 7o} and 
similarly for Qr{xo, to)- Accordingly we define the parabolic distance between sets as 


distpar(A, B) := inf distpar((a;, 7), (y, s)) 
(x,t)GA 
(y,s)eB 


for A,Bc similarly for distpar,G(^i B). 

At a certain point it will be useful to split R” = R”“^ x R. We agree here that we shall 
write a point x G R" as (x', x„) S R”“^ x R; moreover, with B{{x'o) we shall denote the 
ball of R"”^ with radius r and center Xg G R"“^. 
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With S C being a measurable set, xb denotes its characteristic function. If further¬ 
more B has positive and finite measure and / : S —?► is a measurable map, we shall 

denote by 

(/)b = ^ fiy) ^y 

the integral average of / over B. If is a cylinder, B := K x T C then we shall 

denote the slicewise average by 


{f)K{T) := i f{y,T)dy. 
J K 


for almost every r G F. By sup we shall mean possibly the essential supremum, and 
similarly for inf. We shall also as usual denote 

\f{x) - f{y)\ 


osc f := sup f — inf f, 
B ■’ B B 


[f]co.-y{B) ■= sup 


x^y 




Dif := djjdxi, for i G n}, will stand for the partial derivative of / in the Ci 

direction, and Dj will denote f /dxidxj. Here is the i-th element of the standard 
orthonormal basis of K”. By 2* we shall denote the Sobolev conjugate exponent of 2, with 
the agreement that in the case n = 2 we fix the value of 2* as 4, i.e., 

f 2n 

2 * ; n - 2 (2.2) 

[4 n = 2. 

With s being a real number, we shall denote s+ := max{s, 0} and s_ := max{—s, 0}. For 
a vector ^ = (^ 1 ,..., G R", diag ^ denotes the diagonal matrix Finally, 

]R_i_ := [0, 00 ), N is the set {1, 2,... } and No = N U {0}. 

By “equation structurally similar to (1.8)i” we mean an equation of the type dtu — 
div A{Du) = 0 with A satisfying assumptions (1.9) with v, L and g replaced by F, L and 
g. Both will depend on data, while g will satisfy (1.2) and (1.5) with go, go, q 
depending on data and q. 

2.2. Properties of g. Without loss of generality we assume that 


/ 9ip)dp = l. 

Jo 


(2.3) 


Since (1.2) implies that the map r 1 —>■ g{r)r is increasing, while r 1 —>■ g{r)r 

turns out to be decreasing, we have 

minjo;®””^, < g(ar) < max{a®““^, }y(r) 

for every r,a>0; clearly p(0) = 0 and limr_>.oo g(r) = 00 . Since moreover g is strictly 
increasing, it has a strictly increasing inverse function g~^ G C'^(]R+) with 


1 


{9 ^Yir) = -ff-zi, NN 
9 [9 NO) 


for every r > 0 . 


Using (1.2) we then see that also g ^ satisfies an Orlicz-type condition 

for every r > 0 . 


1 


^ {rg N'(r) ^ 


(2.4) 


yi - 1 g N^) flo - 1 
Therefore, anything derived from (1.2) for g holds for g~^ with yo ~ 1 and yi — 1 replaced 
by l/(si ~ 1 ) and l/(yo ~ respectively. 

Define the function G : M+ —>■ M+ as 


G{r) := [ g{p) dp. 
JQ 


(2.5) 
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Clearly G'{r) = g{r) > 0 and G"{r) = g'{r) > 0 implying that G is both strictly 
increasing and strictly convex in (0, oo). Moreover, G(0) = 0 and G(l) = 1 due to 
(2.3). We also define l/G(l/s) = l/G“^(l/s) = 0 for s = 0. It is simple to check by 
integrating the function r i-A rg{r) by parts and using (1.2) that also 


go < 


G'{r)r 


< 9i 


( 2 . 6 ) 


holds true for r > 0. 

Define the Young complement of G as 

G(r) = sup (rs — G(s)) or G{r) := f g~^{p)dp; 
s>0 Jo 


(2.7) 


in our setting these definitions are equivalent, see [33]. Note that the Young’s inequality 


sr < G(s) + G(r-) 


( 2 . 8 ) 


holds true for every r, s > 0 and by (2.4) and the second definition in (2.7) also G satisfies 
an Orlicz-type condition 

gi ^ G^{r)r ^ go 
gi - 1 “ G(r) “ go - 1 ’ 

Now starting from (2.6) and (2.9), we deduce precisely as for g the inequalities 
min{a®“,a®^}G(r) < G{ar) < max{a®“,a®^}G(r), 


(2.9) 

( 2 . 10 ) 


and 

{ 31 _90_ 'I ~ ~ ( 31 3n 'I ~ 

q;3i-'^ , a3o-i >G(r) < G(ar) < max< , aso-i !-G(r) 

for every a,r > 0. These, together with Young’s inequality (2.8), imply for 0 < e < 1 
sr < G(e w s) + G{e~'^r) < eG(s) + c{go, e)G{r). 


Another useful property is 


G 



for every r > 0, 


see again [33] for the easy proof. 

From the second assumption of (1.9) we easily derive an upper bound for A. Indeed, 
when ^ G K" \ {0} we have 


1^(01 < lei f 177^(^01 ds < m f ds 

Jo Jo sKl 

/•l?l Gfifl) 

<c{L,go)J g'{r)dr < c{go,gi,L)—^-, ( 2 . 11 ) 

this holds also for = 0 by our conventions, since .4(0) = 0. Similarly, the first assump¬ 
tion of (1.9) yields 

/•i r\i\ 

{A0,0= {DAis^)^,Ods>c{gi,v)\^\ g'{r)dr > c{go,gi,v)G{\^\). 

Jo Jo 

( 2 . 12 ) 


We define the quantity V], : R" —^ R" by 
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when ^ 0 and set Vg (0) = 0. Clearly Vg is a continuous bijection of R” and, moreover, 
has a continuous inverse by the inverse function theorem. Furthermore, the following 
monotonicity formula holds true: 

(^(a) - -4(6), a - a) > - ap > cif^aa) - ^g(a)P (2.13) 

for a constant c = c{go,gi,i') and for every a, a ^ 1^", see [16, 17]. 


2.3. A concrete example. We give here a nontrivial example of a Lipschitz function g 
satisfying our assumptions - see Remark 1. This example is inspired by [28]. In particular 
we want to demonstrate the possibility that g oscillates between degenerate and singular 
behaviour. Suppose 2nl{n + 2) < po < 5i and set S = {gi — go)/i > 0. Define the 
sequence Sk = 2? for fc G No and the function 


9{s) 


{ sgo-i+<5^ 0 < s < 2 

S 2 fc <a< S 2 k +1 
'S2fc+1 < S < S2k+2 


Clearly g is Lipschitz and it satisfies (1.2). Moreover, (2.3) holds after scaling by a suitable 
normalization constant. We observe that 


lim sup 

S—>-oo 


g(s) 

s 


gi > 2 + (5 

(iff 

go + 2pi > 6) 

gi = 2 + (5 

(iff 

go + 2pi = 6), 

gi < 2 + (5 

(iff 

go + 2pi < 6) 


lim inf ^ 


s—>-oo s 


go > 2 - (5 

(iff 

2go + gi > 6) 

O 

II 

to 

1 

(iff 

2go + gi = 6) 

go < 2 - (5 

(iff 

2go + gi < 6) 


By taking po = 2 — ^= 2 + ^ we obtain a particularly interesting case, that is, we 
have liminfs_>.oo g{s)/s = 0 but limsup^_^^ 9 {s)/s = oo. Furthermore, if we consider 
the function 

' im- 

we find similar behaviour as s 0. This is to say, we can build a structure function g (and 
accordingly a vector field ^ as in (1.1)) that, forG N, along the sequence the 

function g[s )/s is at the same time as large and as close to zero as we wish, and therefore 
it does not enjoy any monotonicity properties. This gives a clue about the difficulty of 
the application of De Giorgi-type methods, in particular when they have to be matched 
with intrinsic geometries: note that the expressions of the type G{s)/~ g{s)/s appear 
already in the energy estimate for (1.1), see Lemma 2.3. On the other hand, when the 
quantity g{\Du\)/\Du\ is known to be under control, then the equation becomes treatable, 
see for instance Proposition 3.4 and in particular (3.8). 


2.4. Orlicz spaces. For G as in (2.5), a measurable function it : A —R, A C R^, k gN 
belongs to the Orlicz space L'^{A) if it satisfies 


G(|it|) dx < oo. 


The space L^{A) is a vector space, since G satisfies the A2-condition (2.10), and it can be 
shown to be a Banach space if endowed with the Luxemburg norm 


:= inf 
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A function u belongs to L'^{A), if u G L'^{A') for every A' (s A. If also the weak 
gradient of u belongs to L^{A), we say that u G W^’^(A). The corresponding space with 
zero boundary values, denoted Wq’‘^ (A), is the completion of (7“ (A) under the norm 

ll'wllivi.G(^) := ||u||iG(^) + \\Du\\lgi^a)- 

We denote by V^{Q,t) the space of functions u € L^{Q,t) H for 

which also the weak spatial gradient Du belongs to The space is also 

a Banach space with the norm 

l|u||vG(o,j,) := ||u||iG(o,j,) + \\Du\\lg(^q^). 

Moreover, we denote by V'Q^(flT) the space of functions u G V^(flT) that belong to 
Wq’ (n) for almost every t G (0,T), while the localized version is defined, as 

above, in the customary way. We also shorten 

:= L°° (0, T; L‘^{n)) n V^inr) 

and similarly for the localized and the zero trace versions. We shall moreover denote 
V'^’P{nT), forp > 1, the space for the choice G{s) = s^. 

2.5. The concept of solution and consequences. We fix here the notions of solution em¬ 
ployed in this paper. 

Definition 1. A function u is a weak solution to (1.8)i in a cylindrical domain /C C 
with the vector field A satisfying the assumptions (1.9), if u G (^) it satisfies the 
weak formulation 

I [ — udtT] + {A{Du), Drj)] dx dt = 0 (2.14) 

JK 

for every test function rj G {1C). If instead of equality we have the <(>) sign for every 
nonnegative 77 G C^{IC), we say that u is a weak subsolution (supersolution) in 1C. 

Definition 2. A function rt is a solution to the Cauchy-Dirichlet problem (1.8) if u G 
c°(nT) is a weak solution to (1.8) 1 in CIt and moreover u = pointwise on dpCtr- 

A very useful formulation, equivalent to (2.14), is the one involving Steklov averages. 
Indeed, the mild regularity of a solution does not allow us to use it as a test function. 
Furthermore, it is sometimes useful to have a weak formulation allowing for test functions 
independent of time, or test functions possibly vanishing only on the parabolic boundary 
of a cylinder. Apart from mollification, the possible way to have such properties involve 
the so-called Steklov averaging regularization of a function: for f : 1C = V x {ti, 12 ) ^ ^ 
measurable and 0 < \h\ ^ 1 appropriate, it is defined as 

1 /■* 

fh{x,t)-=T f{x,s)ds for {x,t) GV X {ti + hA 2 )] 

h Jt-h 

note that we employ the backward regularization. If f G for some q > 1, then 

//j —>■ / in L'*(T>x ^2)) forevery e > 0; the same holds in the spaces. Moreover, 

if / G C°{ti,t 2 ', L^{TC>)) then fh{-,x) /(u x) in L^{T>) for a.e. t G {ti + e, ^2) and for 

every e > 0. 

At this point it is quite easy to infer the following slicewise formulation for weak solu¬ 
tions (see [14]) using density arguments with respect to the spatial variable: 

[ [dtUhi-,T)ri + {[A{Du)]h{-,T),Dr]}] dx = 0 (2.15) 

Jv 

for every 77 G VFq’*^( 22), almost every t G {ti + h, ^2), and h > 0 such that the functions 
are well defined. Similar results hold also for weak super- and subsolutions. 
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Proposition 2.1. (Comparison principle) Lef /C ■='Dx (<i,< 2 ) C flrandletu € C^{Kf) 
be a weak subsolution to (1.8)i and v C C^(K^) a weak supersolution to (1.8)i in K,. If 
u < V on dpK,, then u < v in K. . 

Proof. For £ > 0 fixed define (pe{t) ■= (^2 — £ — 0+ and test (2.14) formally with 

V = {uh - Vh - s)+P6- 

Note that p is compactly supported in /C due to the continuity of u and v and the fact that 
u < u on dpQ. Subtracting the Steklov version of the variational inequality of v from that 
of u and integrating over (^1,^2) yields 

I dtiuh - Vh)r] dx dt + j {[A{Du)]h - [A{Dv)]h, Dp) dx dt < 0. 

JK JK 

By the monotonicity of A, Lemma 2.13, we have 

[ {[A{Du)]h - [A{Dv)]h,Dp) dxdt 

JK 

—>■ / {A{Du) — A{Dv), Du — Dv))ipe dx dt > 0, 

J Kr\{u>v+e} 

and for the parabolic term we obtain using integration by parts 

I dt{uh - Vh)(uh -Vh- e)+Lpedxdt = [ {uh - vt - e)\dt(pedxdt 

JK ^ JK 


K 
t2 — e 


uriA 


— V — £)^ dx dt 


as h —>■ 0. Combining these gives 

pt2 — e 


ri2—s r 

/ / {u — v — e)\dxdt<Q, 

Jti Jv 


which implies u < v + e almost everywhere in 2? x {ti,t 2 — e)- Since this holds for every 
e > 0 and u,v G C{Q), the result follows. □ 

Observe that the uniqueness of a solution to the Cauchy-Dirichlet problem (1.8) follows 
immediately from the previous result. Moreover, we have the following corollary. 

Corollary 2.2. (Maximum principle) Let K. C fir and let u G C{Kf) be a weak solution 
to ( 1.8)1 in K.. Then 


in KA and, moreover. 


inf u < u < sup u 
SpK- d„K 


sup |u| = sup |m|. 
K dpK 


We recall the following standard energy inequality for local weak solutions. We give it 
in a more general form for future reference. 

Lemma 2.3 (Caccioppoli’s inequality). Let K. :='D x (^1,^2) <£ and let u be a weak 
solution to (1.8) in K,. Then there exists a constant c = c{go,gi, v, L) such that 


sup f [(u — (•, r) da; + f G[\D{u — k)±\)ip^^ dx dt 

r&{ti,t2)JV JK 


< 


[{u-k)\ip^^]{-,ti)dx 


G{\Dtp\{u - k)±) +{u- k)\.\dtip\ 


dx dt 


Jv JK ^ 

for any k gM. and for every ip G vanishing in a neighborhood of dD X (fi, t 2 ) 

and with 0 < tp < 1. The same inequality but only with the “+” sign holds for weak 
subsolutions. 
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Proof. Fix tf G as in the statement of the Lemma, call w := ±(u — k)± 

and choose 77 = as the test function in (2.15). Then we integrate over (fi,T) for 

T G (ti,t 2 ) to obtain 

/ dtUhWh(p^^X{ti,r)dxdt + {[A{Du)]h,D{wh(p^^))x{ti,r)dxdt = 0. (2.16) 

Jk Jk 

Integration by parts gives 

[ dtUhWh^^^X{ti,T)dxdt = ^ f f dt{wl)ip^^ dx dt (2.17) 

Jk. ^ Jti Jv 




dx 


^ It J^^hdtiv^")dxdt 




1 dx 

2 Jv 

as h —?► 0. For the elliptic part we have by (2.12) 

f {[A{Du)]h,D{whip^^))x{tuT)dxdt 




dxdt 


' J'D 


{^A{Du), Dw) dx dt + gi f f (^A{Du), Dip) w ^dxdt 

Jt\ J 7 ? 

> Cl f f G{\Dw\)p^^ dx dt — gi f {^A{Du), Dp'j w dx dt 

Jti Jv Jk 

where ci depends on go,gi, v. Furthermore, by (2. 1 1), Young’s inequality with e G (0,1) 
to be chosen and the properties of g we obtain 


gi / (A(Du) , Dp) w ^dxdt 
Jk 

< ec2 [ G 


<gi |.4(Dr(;)||i7(^||?u|(^®^ ^dxdt 
Jk 

(dx dt + c{e) [ G{\Dp\\w\) dx dt 

JK \ \dJw\ ) 

< ec 2 f G{\Dw\)p^^ dx dt + c{e) ( G[\Dp\\w\)dxdt, (2.18) 

Jk Jk 

where C 2 depends on gQ,gi,L and c(e:) depends on poj 5 ii ^ as well as on e. Now, com¬ 
bining (2.17)-(2.18) with (2.16) yields 

^ f vAp^^ dx ~^ [ w'^dtip^^) dxdt + Cl f f G{\Dw\)p^^ dxdt 
2 Jv t=ti 2 Jk Jti Jv 

< ec 2 I G{\Dw\)p^^ dx dt + c{e) f G{\Dp\\w\)dxdt. 

Jk Jk 

We conclude by taking the essential supremum with respect to r G (ti,t 2 ), choosing 
e G (0,1) such that ec 2 < ci/2, reabsorbing the term on the right-hand side and recalling 
the definition of w. 

The proof for subsolutions is very similar, taking into account that the test function 77 
must be nonnegative. □ 

2.6. The geometry of the problem. In order to understand the equation, the first thing 
we want to stress is its scaling. Suppose u solves the model equation (1.1) in Qi = 
Bi X (—1,0) and let k > 0. Then the function 

'X — Xq 1 K' 


i{x,t) :=Ku(^-^,^G(^)(f-fo)) 


Q'^{xo,to) ■= Br{xo) X (to - ’^ 0 ) 


solves in 
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the equation 


where 


/ g{\Du\) \ 


\Di 


Du = 0, 




(2.19) 


( 2 . 20 ) 


The function g has the same structure as g, in the sense that it satisfies (1.2) exactly with 
parameters go and gi and moreover, we have G(l) = 1, where 


G(s) := ^ g{a) da = [g(^)] g(^s) . 


Conversely, if we have a solution w to (1.8) in QJt, then 


'j{x, t) := —w( 
K V 


xo + rx, to + K 


G 


(7) 


1 -1 


t 


solves (2.19) in Qi with g as in (2.20). In case we consider the general equation (1.8), the 
same scaling argument holds if we consider the vector field 




which satisfies the structural conditions (1.9) with g replaced by the function g. 


2.7. Other auxiliary results. The following Lemma encodes the self-improving property 
of reverse Holder inequalities. We take the form proposed in [27, Lemma 5.1] with slight 
changes in order to meet our purposes. 


Lemma 2.4. Let g be a nonnegative Borel measure with finite total mass. Moreover, let 
7 > 1 and {crQ}o<CT<i bs a family of open sets with the property 

<y'Q C crQ C IQ = Q 


whenever 0 < cr' < cr < 1. Ifw € Lf{Q) is a nonnegative function satisfying 


w 


27 


\ 1 /( 27 ) 


G'Q 


dpj 


< 


Co 


a — a' 


'(tQ 


uP' dp 


1/2 


for all 1/2 < a' < a < 1, then for any 0 < q < 2 there is a positive constant c = 
c(co, 7, q) such that 


,,27 


1 /( 27 ) 


’(jQ 


dp 


< 


(1 - cr)« 


for all Q < a < 1, where ^ := 

The next one is a classic iteration Lemma. 



Lemma 2.5. Let f : [i?, 2R] — >■ [0, c») be a function such that 

1 A 

(/{t) < -4>(s) + 7-rj + B for every R < r < s < 2R, 

2 [s — rjP 

where A,B >1 and /3 > 0. Then 
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3. A PRIORI Lipschitz estimates 


In this section we impose on u an additional regularity assumption and prove intrinsic 
estimates for the gradient of u. To be precise, we shall suppose 

u^DuGC^nr), uGLUO,T;Wt\n)). (3.1) 

This is to say, we shall prove the estimates of this section as a priori estimates, leaving 
to Section 4 the approximation procedure which will explain how to deduce the desired 
estimates without the additional assumption (3.1). Notice that the continuity of u and Du 
allows us to treat their pointwise values. Due to the assumed extra regularity it will be 
possible to differentiate the equation; this will be done by showing that the function 

V := \Du\‘^ (3.2) 

is a subsolution to a similar equation. 

Lemma 3.1. Let u be a weak solution to (1.8)i in LIt and, moreover, assume that the 
regularity assumptions (3.1) hold. Then v is a weak subsolution to 

dtv — div{DA{Du)Dv) = 0 in LLt- (3.3) 


Proof. Formally, the idea is to differentiate equation (1.8)i with respect to Xj for j = 
1 ,... , n , then multiply by DjU , and finally sum over j. To this end, let 0 < (^ e C-COt), 
and test (2.14) with 

r ] = — Dj { DjU ( p ). 

This choice can be justified by using Steklov averages, as done previously in the paper; we 
shall proceed formally. Integration by parts yields 


0 = - 


/ udt{—Dj{DjUip))dxdt+ / {A{Du), D{—Dj{DjU(p))) dx dt 

J O-f D Cl-J’ 

/ dt{Dju)DjUip dx dt + / {DjA{Du), D{DjUip)) dx dt 

— i j \Dju\^dt(pdx dt + ^ I {DA{Du)D(^\Dju\^),Dip) dxdt 

2 Jnj. 2 

+ / {DA{Du)DDjU, DDju)(p dx dt. 

J 


Now, since 


[ {DA{Du)DDjU, DDju)ip dx dt > V f j^ DDjuf (p dx dt > 0 

J J vl'i' 

by (1-9) summing up over j = 1,..., n leads to 

— / \Du\^dtp dx dt / {DA{Du)D\Du\'^,Dp) dxdt <Q. 

J J 

This proves the claim. 


□ 


Next we prove a Caccioppoli inequality of porous medium type for the function v. 

Lemma 3.2. Let u be a weak solution of (1.8) in Llx and assume that (3.1) holds. Let 
K. :=D X (fi, t 2 ) <£ LIt and fc G R. Then there exists a constant c = c(y, L) such that 

sup [ [{v - k)\p'^]{-, t) dx + [ - k) + \'^p'^ dxdt 

r£iti,t2)Jv Jk 

< cj^iv - k)l + \dtP^ dxdt 

for every p G C°°{IC) vanishing in a neighborhood of dpIC. 
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Proof. We can take 

r]={v- fc) + (/7^X(ti,r) 

for r € (^1,^2) as the test function in the weak formulation of (3.3), up to a regularization 
similar to the previous ones. For the parabolic part we have 

— f f V dt{{v — dx dt = ^ j f dt{v — k)\ip^ dx dt 

Jti Jv 2 

= \j [(^-j {v - k)ldt>p'^ dxdt. 

The elliptic term can be estimated from below by using the assumptions ( 1 .9) and Young’s 
inequality with e = i/l{2L). This gives 

f f {^DA{Du)Dv, D(^{v — k)+(p^)) dx dt 
Jt\ J'D 

= f f {DA{Du)D{v — k)+, D{v — k)+)(p^ dx dt 
Jti Jv 

+ 2 / f {DA{Du)D{v — k)+, Dip) {v — k)+(p dx dt 
Jti Jv 

^ r r g{\Du \), 

Jti Jv \Du\ 


> V 


-\D{v — /c)+| dxdt 


- 2L 


9i\Du\) 


u 

~ 2 


ti Jv \Du 

9i\Du\) 


\D{v — k)+\\Dp\ {v — k)+pdxdt 


11 J TJ) 


\Du\ 


\D{v — fc)+| dxdt 


-civ,L) f f {v - k)ldxdt, 

Jti Jv \Du\ 

Jjiv - r) dx + uj^^ dADu \), 


and thus, we obtain 


-\D{v — k)+\ p^ dx dt 


< c 


' ti 'J TJ) 

9{\Du\) 

\Du\ 


\Dp\‘^ {v — k)'^ dx dt + / {v — k)'^\dtp\ dx dt. 

Jk 


Since r G (^1,^2) was arbitrary, the result follows. 


□ 


Combining the previous lemma with Sobolev’s inequality leads to the following esti¬ 
mate. 

Lemma 3.3. Let the assumptions of Lemma 3.2 be in force. Then there exists a constant 
c = c(n, pi, v, L) such that 


i 


- k)^^p'^'^ dx dt 

\Du\ 


|T>p/”(f2 - fi)^ - k)l(^jJj^^\Dp\'^+ \dtp\^ dxdt^ , (3.4) 


<c| 

where - recall (2.2) - 

2 

7 := 2 - — > 1. 
2 * 

Proof. By Holder’s and Sobolev’s inequalities we have 
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< 


^2 — 
1 


1 f 

~Ju It 


I 'D 
t2 


h - ti Jt^ 


_ k)i^mv - dxdt 

2/2* 


£(( 


\Du\ 

g{\Du\)\U-2 


\D% 


sl/2 Y 

j (ti — k)+ip j dx 1 

-f (U- fc) + /9^ 

Jt> 


dx 


1 - 2 / 2 * 


dt 


1 - 2 / 2 * 


< c{n)\'D\^^" f sup I [{v — k)'^ip'^]{-,T) dx 
\T&{ti,t2)Jv 


dxdt. (3.5) 


A straightforward calculation yields 


D 


isilDuDY// ,, 5 


(ti- fc)+ / \Du\g'{\Du\) ' 

4u 'v g{\Du\) 


(W) 


and thus, integrating and estimating the first term using Lemma 3.2 yields 

'g{\Du\)\h '' ^ 




dx dt 

< c 


i {v-k) 
JK 


( 9{\Dn\) 

\ \Du\ 


\Dlp\'^ + \dtg)\ I dxdt 


where the constant c depends only on gi^iy, and L. From Lemma 3.2 it also follows that 

I \(v - 


sup f [{v - k)_^_ip ]{-,T)dx 
Jt) 

< ciU - I 


{h - h) - k)\ l-P/gp + dxdt] 


therefore, by inserting the previous two inequalities into (3.5) we obtain (3.4). 


□ 


Next the aim is to prove an intrinsic reverse Holder’s inequality. To this end, let 
Qp{x[), to) C Ht, let a > 1 be such that 

A > - sup \Du\, (3.6) 

4 Qp(xo.to) 

and set 

0 ■= ^ 

^ ■ A ■ 

We introduce the intrinsic cylinder 

Qp = Qpixo,to) :=m.m{l,exy^^Bp{xQ) X (fo - min{l, fo)- (3.7) 
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Note that we have the alternative expression 

'Bp{xo) X {to-9x^p'^,to), 0x>l 

0l^^Bp{xo) X (to - O<0x<l, 


Q.= 


from which we easily see the analogy with the intrinsic geometry used to handle the par¬ 
abolic p-Laplacian, recalling that in this case g{s)/s = 5^“^ and A is “dimensionally 
comparable” to \Du\. Observe that we clearly have Qp{xQ, tg) C Qp{xQ, tg) in any case. 

Lemma 3.4. Let u be a weak solution to (1.8)i in LIt, assume that (3.1) and (3.6) hold 
and let q > 0. Then there exists a constant c = c{n, pi, i/, L, q) such that 


( -/ {v — k)'^ dxdt] < c 


\ 1/(27) 


/ (v-k) 

JoL 


— dx dt j 


\ 1/9 


for every k > 

Proof. Let 1/2 < cr' < cr < 1 and choose a cut-off function ip G C°°(aQp) vanishing in 
the neighborhood of dp(aQp) such that 0 <(/j< 1 , (/J=lin (j'Q\ and 


\Dt\ < 


p(cr - cr') 


{1,0 a} 


- 1/2 


IdM < 


{1,0A-^}-^ 


— cr')2 


Observe that by the inclusion Qp(xo, to) C Qp{xo, to) and (3.6) we have 

\Du\ < 4A in Qp. 

Moreover, we have \Du\ > A in the support of {v — k)+, since k > X? and v = \Du\^. 
Thus, by using the properties of g we obtain 


Ici / 9{\Du\) ^ 
4^^ - -\D^ - 


(3.8) 


in QA n {t; > k}. Now Lemma 3.3 yields 


-f {v — k)^ dx dt < c{n) 0 -f ^^j——^{v — k)'^(p‘^'^dxdt 

I 2jn / \ 7 —1 

< min{l,6'A}^/^.B^p min{l, ff^/^Kcrp)^ j 


< ^ min{l, 0 a} iiiin{l, 0 A ^ 

'0x min{l,0A}“^ + min{l,0)/^}“ 
p2(cr — cr')2 


■ / {v-k) 


^ dx dt 


(cr - cr'; 


This is to say 


/ r A 1/(27) r f r 


where the constant c depends only on n, gi,v, L. 
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Next we use Lemma 2.4 with w = {v — k)+ and dfi = dx dt. This gives for every 
0 < g < 2 a constant c = c{n, gi, u, L, q) such that 

/ r \ 1/(27) / p \ 1/9 

yj- {v — k)'^dxdtj < cyj- {v — k)'^dxdtj ; 

the case q >2 now follows from Holder’s inequality. □ 

Iterating the previous result yields the following pointwise estimate. 

Proposition 3.5. Let u be a weak solution to (1.8) in Llx and assume that (3.1) holds. 
Then for every g > 0 there exists a constant c = c{n, gi, v, L, g) such that 

/ p \ 1/(29) 

\Du{xQ,tQ)\ < A + c ( -f {\Du\^ — dxdt j 

Q^{xo,to) ) 

holds for every X satisfying (3.6). 

Proof. The idea is to apply De Giorgi’s iteration method with the aid of Lemma 3.4. Let 
us first consider the case 0 < g < 2. To this end, choose for j G No 

Pj = 2-1 p, kj = A^ + (1 - 2-i)d, 

where d > 0 is to be determined later. Observe that po = p, ko = X^, and pj decreases 
to zero and kj increases to A^ + d as j tends to infinity; clearly kj > X^. Denote Qj := 
QpAxo,to) and 


1/9 


(v — dx dt I for j G Nq. 


1 /( 27 ) 


’’ (!<’ 

By Lemma 3.4 we have 

(-I- {v — kj)'y dx dt 

and since kjj^i > kj implies 

(u - kjfy > (%+i - kjf"<-\v - kj+i)lx{v>k,+r}, 


< c 


J" A ~ %)+ dx df^ 


1/9 


we obtain 

Yj+i < 


f p \ iil<i 

■\^j {v- kjfA dx dtj < c*d-^2^iY^+^ 


{kj+i - kjY 

for every j G Ng, where j3 := 2y/q — 1 > 0 and c* = c*(n, pi, u, L, g). Then a standard 
hyper-geometric iteration lemma implies Yj —0 as j — 00 , provided that 

Yo < (2c*)"Pd 

and this can be guaranteed by choosing 

1/9 

^ ,dcr dt 


d = 


{2c*y(-f 


’ Qpixo . to ) 

Now Lebesgue’s differentiation theorem yields 


A — A^)^ dxdt j 


(u(a:o,fo) — (A^ + d)) = \im (-f (v—(X^ + d)Ydxdt\ < lim Y, = 0, 

+ J-i-OO \Jq, + ) 

which implies, recalling the choice of d, 

\ 1/9 


v{xo,to) < X"^ + c(^ [v — X^Y dxdt^ 

\JQYxo,to) J 


The case g > 2 follows again by Holder’s inequality. 


□ 
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4. Approximation 

In this section we regularize the equation in order to apply the results of the previous 
section and show that the gradient of the solution to the regularized equation is uniformly 
bounded. Then all we have left to prove is that the approximating solutions converge to a 
function that solves the original equation. 

To this end, define for e e (0,1) 

:= (</>. * Am + e(l + (4.1) 

where ipeiO = ^ is a standard mollifier with cjidx = 1. That is, we 

mollify the vector field A and perturb it with the nondegenerate gi-Laplacian, where 'gi > 
max{( 7 i, 2}; we can take for example gi := gi + 1. It is straightforward to see that Ae 
satisfies (1.9) with g replaced by 

^g(s) := + e(l + s)9i-2s ( 4 . 2 ) 

s + e 

and L, v replaced by L = c(n, gi)L, v = vjciri^ gi), see also Paragraph 6.1. Now the key 
point is that Og^ can be bounded independently ofe. Indeed, we have 

flO - 1 < Ogm < 5l - 1. 

where po := minlsoi 2}. Note that g^ also satisfies the lower bound in (1.5), since ge{s) > 
g{s)/2 for s > 1 . 

Let Ug G 1/2,91 (fly) n be the solution to the Cauchy-Dirichlet problem 

{ dtUe — diy AeiDug) =0 in Ot, 

Ue = Ip on dpflr', 

for existence and uniqueness of such solutions see for instance [24]. Since 

e(l + s)9i-2 < < £M (1 + s) 9 i- 2 ^ 

s e 

in addition to satisfying ^g-ellipticity and -growth conditions analogous to (1.9), the vector 
field Ae also enjoys nondegenerate p-Lap\acian growth conditions with p = gi. Hence, by 
standard theory, Ug satisfies the assumption (3.1), see [14, 27]; therefore the results of the 
previous section are at our disposal for u = Ug- Note that all the constants will turn out to 
be effectively independent of e. 

Let us then show how to apply the result of the previous section in order to locally 
bound the gradient of the approximating solution uniformly in terms of e. Here we also 
prove an estimate that, once convergence is established, leads to (1.10). Observe that the 
assumption (1.5) is crucial in this proof. We shall shorten ||^||yoo = ||V'||L“(apOT)- 

Proposition 4.1. Let Ug be a solution to (4.3) and let K. (s LIt- Then ||i9M£||yoc,(^) 
is bounded by a constant depending on data,e,Ci, and distpar(Op^T,^), but 

independent of e. 

Proof. Let us consider a standard parabolic cylinder Qar = Qar{x*, t*) C LIt and a 
subcylinder Qp{xQ,tf) C Q 2 R- Moreover, let A > 1 be such that 

A>i sup \ Dug \. (4.4) 

^ Qp{xa,to) 

We divide the proof into two cases depending on which term of gg dominates at A. 
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Case I. Assume 


Setting 


ei := 


5(A±£) + 

A+e “ ^ ' 




A 


A “h £ 


we clearly have 

e(l + A)si -2 < 0 ^ < 2e{l + A)si- 2 _ 
By applying Proposition 3.5 to with q = pi/2 we obtain 


\Dueixo,to)\ < X + c( -f (l-Duep - da: dA 

\jQ^(xo,to) J 

max{l, 6 >|} f 

Jc 


^ A -|- c 




dxdt 


1/si 


(4.5) 


since Qp(a:o,to) C Qpixo,to). 

We further distinguish two cases; in the case when > 1 we get 


max{l, 9\} 
min{l, d |}"/2 


91 < 2£(1 +A)3l-^ 


while when 0 < < 1 we have 


max{l, 9f} 
min{l, 9^}"^/^ 


(0e)-n/2 < = £:(ei+2/"(l + 


-n/2 


(4.6) 


in both cases we have used (4.5). Since 


g(^ + g) 

- (i + A)3i-2(A + e) 

> c^(A + e)^+'”^(l + A)"”®^ 

> q (1 + A) 2 - 9 i+“”{- 4 /(n+ 2 ). 0 } ^ 


by (1.5) and the fact that A > 1, plugging this estimate into (4.6) yields 


a direct computation shows indeed the relation between the exponents. Hence we have 

In f + \ I tl I \\l-min{e(n+2)/(2gi),2/gi} 

|iA-«e(Xo, fo)| < A + C(1 + A) 



by Young’s inequality; we also used gi > 2. 
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Case II. Suppose then that 


Here we have 


A + £ ^ ’ 


A A “h £■ 


and again by Proposition 3.5 

\DUe{xo,tQ)\ < A + c[ / {\DUe\'^ - X^y dxd^ 

When 9\> 1, choosing q = 1 leads to 


\DUe{xo,to)\ < A + c 
^ A “t" c 


/ g(A + e) 
^ A + £ 

/ 

J Qp{xo,to) 


■/ {\Due\^ - Xy^dxdtj 

jQp(xo,ta) J 


g(l-P^el) 

\DUe\ 


(IZJueI^ — A^) dx dt 


< 


X + c(-h G{\Due\)x{\Du,\>i} dx dA . 

\jQp{xo,to) / 


The second inequality stems from the fact that 

1 


— (A + £) ^ l-Drtgl < 4(A + e) 


(4.7) 


in the set Qp{xo,to) fl {|77Me| > A} by (4.4), while for the last one we used (2.6) and the 
fact that A > 1. 

In the case 0 < < 1 we choose q = e(n + 2)/4 and use (1.5) and again (4.7) to 

obtain 

|£)Ue(a;o,fo)| < A + c^^^y^^^ J ^ {iDuel"^ - X^Y^ dx dt^ 


< 


Qp{xo,to) 

X + c(-f 

\J Q„(xn,tn) 


-e)f+2(j 


— A -|- c 


/ 

Jc 


Qp(xo,to) 


\DUe\ ^-+^^^X{\Due\>l}d^dt 


\ e(n+2) 

I 


note that 


(- 


< 


n — 2 \ n 


n + 2 

Therefore in both cases we have 


A + c(+ G{\Due\)x{\Dup\>i}dxdt\ 

\jQp{xo,to) / 

\n /4 \n n+2 


«(". + 2) 


.n + 2 


)? 


, n — 2 

= 14-+e. 


|Dne(a;o,fo)| < A + c( 4 G{\Due\)x{\Dup\>i} dx dtj 

\jQpixo,to) J 

Combining Cases I and II and denoting fj := max {1, } yields 

|£>UE(xo,fo)| 


n + 2 
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<2X + c( -f (G{\Due\)x{\Du,\>i}+£\Due\^^\dxdt] +1 

V J Qp(xo.io) V J J 

{n+2)fi / f \ ^ 

<2A + c|^—j Ge{\Due\)dxdtj +1, (4.8) 

since 

G{s) < —g{s + e)(s + e) < < —55^(5) < lliG'e(s) 

9o 9o s + e go go 

for s > 1 and trivially 

< e(l + < giGe{s). 

The constant c in (4.8) depends only on data, e, ce. 


Let us now choose two intermediate cylinders Qn C Qr ^ Qs C Q 2 R and fix 
1 S — T 

^ 1 + ^II^We||L=o(Q^) < 00, {xo,to)&Qr, p := ^—> 0. 

Clearly Qp{xo, to) C Qs so that (4.4) holds. Then (4.8) implies 


\\DUe\\L’^(^Q^) < 2^\DUe\\L°°{Q^) 

/ R \ (n+2)ij / f \ ^ 

+ C^-j Ge{\DUe\) dx dtj +3. 

Now, by choosing (j){r) = \\Due\\L°°{Qr)’ iteration Lemma 2.5 gives 

\\Due\\L'^(Qj,) < c( ^ [Ge{\Due\)+ i\dxdt\ . (4.9) 

\JQ2R / 

At this point, in order to get rid of the dependence on e on the right-hand side, the idea is to 
use the Caccioppoli inequality of Lemma 2.3 to translate the dependence on Du^ to one on 
Me, and the latter in turn into a dependence on ijj- Indeed, take (f € G°° {Qar) vanishing in 
a neighborhood of dpQiR such that 0<ip<l, (p = lin Q 2 R, and \Dtf\’^ + \dtif \ < c/R?. 
Since 

SUp|Ue| < sup |Me| = SUp |^/>| < ||V’|U“ 

Qar dpQ.T dpO^T 

by the maximum principle. Corollary 2.2, we can estimate by Lemma 2.3 



Ge{\DUe\) dx dt < 


c-f [Ge(|i9(p||Me|) 
^ Qar 


+ Mel^t'/’l] dxdt 


< c 





= c(data,e,Q, ||V’|U«=,i?)- 


(4.10) 


Note that the constant does not depend on e. Therefore we conclude the proof of the 
Proposition, modulo a standard covering argument. □ 


4.1. A uniform interior modulus of continuity via Lipschitz regularity. In this section 
we prove that the approximating solutions are equicontinuous in the interior of the 
domain; in particular we shall show their equi-Lipschitz regularity with respect to the 
parabolic metric. 

Proposition 4.2. Let be a solution to (4.3). Then G Lip(l, l/2)(ilT) locally, uni¬ 
formly in e; this is to say, for every subcylinder K, g LIt there exists a constant c depending 
on data, e, ci, ond distpar(9pflT, Ai!) such that 

\Ue{x,t) - Ueiy,s)\ < C distpar ((a:, f) , (?/, s)) 


(4.11) 
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for every (x, t), (y, s) € 1C and for every e € (0,1). 

Proof. Fix an intermediate set Kf such that JC <£ K.' (s CIt and 

distpar(5, = distpar(/C,9pf2T)/2 =: d/2 

for every z G dpKf. Take also a cylinder Qr{xo,to) C Kf with (a;o,fo) C 1C', this will 
happen for instance if r < d/2. Since Du^ is continuous, by applying the divergence 
theorem and using the bound for Ae in (2. 1 1) we infer 


£ 


Brixo) 


Ue{-,T)dx 


=^ / 7 

T = ti ^ Jti J dBr(xo) ' 

^ Jti J d 


^ ' d-H^-^ dt 


F - a;o| 

geiDue) dn^-^ dt 


I dBr{xo) 

for all to — < ti < t 2 < to, where stands for the (n — 1)-dimensional Hausdorff 

measure. We thus estimate 


Te{to — r^,to) 


sup 


4 Ue{',T)dx 

t2 

\ J Br{xa) 

T=ti 


^ c 


/ 

Jdl 


g,{Du,) dn^-^ dt 


^ J tg — r^ J dBr{xo) 

< cr(l + 

<Cr{l+\\DUe\\L°-(K')Y^ 

Now by Proposition 4.1, in particular by (4.9)-(4.10), we have 

,e, ||^||r~,a) r. 


osc iue)Br.{xo){T) < c(data,c^,e, ||V'||L~,d)' 

Te(to—r^,to) 

At this point we simply split for (xi, fi), (a;2, <2) G Qr(xo,to) 


(4.12) 


-{Xi,ti) - Ue{x2,t2)\ < 

Ue{xi,ti) - f 


Ue{',ti)dx 


J Br 

(xo) 


f Ue{',ti)dx - 4 

Ue{-,t 2 ) dx 

+ 

Ue{x2,t2)-f Ue{',t2)dx 

J Br(xo) J 

Br(xo) 


J Br{xo) 


+ 


While in order to bound the second term we shall use (4.12), the first and last terms can be 
estimated using the mean value theorem as follows: 

;{Xi,ti) - + Ue{-,ti)dx < f \ue{Xi,ti) - Ue{x,ti)\dx 

J Br-(xo) J Br(xo) 

< 2r \\DUe\\L’^(K.'), 

for i G {1, 2}. Therefore, using again Proposition 4.1, we have 

osc Ue < cr (4.13) 

Qr{xa,to) 

with c as in (4.12), in particular not depending on e. To conclude the proof, for (xi,fi), 
{x 2 ,t 2 ) G /C, we simply check whether distpar((a:i, fi), (x2, f2)) < c?/4 holds true or not; 
if so, then there exists a cylinder Qr(xo,to) with r = distpar((xi,fi), (x2,t2)) such that 
(xi,fi), (x2,t2) G QrixoAo) and we can apply (4.13) that directly yields (4.11). If on 
the other hand distpar((xi,fi), (x2,f2)) > d/f then, again simply using the maximum 
principle, we have 

I ^ ! M ^ oil II ^ distpar((a:i, fl), (X2, < 2 )) II , II ,A-tA\ 

\Ue{x,t) -Ue{y,s)\ < 2\\Ue\\L’^(a.j.) < 8- - -llV'lk-; (4-14) 

□ 


the proof is concluded. 
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Remark 2. Notice that, tracking the dependence on d of the constant in Proposition 4.2 and 
in turn the dependence on R of estimate (4.10), and also slightly modifying the previous 
proof, we deduce that estimate (4. 1 1) can be rewritten as 

\ ue { x , t ) - Ue{y,s)\ < distpar ((x, f), (y, s)), (4.15) 

^z,w 

for an exponent 7 = 7 ( 72 , pi,e) > 1 and a constant c depending only on data,e, q, 
||iA||l“, with 2 = (a:, t), w = (y, s) and accordingly 

dz,w ■■= min{distpar(2,9pf2T),distpar(w,9pf^T),l}- 

Indeed, if distpar(z, w) < dz.iu/S, then we can apply the argument in the first part of the 
proof of Proposition 4.2 with r = distpar(2:, w) to get (suppose s <t) 

Q 

\Ueiz) -Ue{w)\ < OSC distpar(z, tu), 

Qr{z) dz,W 

where 7 = 51(51 — 1)5, since we have Qr{z) C C fir and so 

\\DUe\\L-^(Q^(z)) < ll^“ellL“(Qd,,„/8W) ^ 

Clz,w 

The case where dz^w < 8 distpar(z, w) can be approached exactly as in (4.14). 

5. Continuity at the boundary 

In this section we prove that the solution to the approximating problem (4.3) is continu¬ 
ous up to the boundary independently of e by building an explicit barrier. We do not want 
to enter the details of the theory and the general relation between existence of barriers and 
regularity of the boundary points; the interested reader can see the nice paper [24] for the 
evolutionary p-Laplacian, while [21, 32] summarize the results in the elliptic setting. 

We shall begin with the proof of the continuity at the lateral boundary; here we shall 
give all the details needed. For the continuity at the initial boundary we shall however only 
sketch the proof, which on the other hand is very similar and easier than the lateral case. 
Again, we will prove the existence of a uniform (in the sense that it will be independent 
of e) modulus of continuity for in the last section we shall show that this modulus is 
easily inherited by the limit of u ^- 

Let us begin with the construction of an explicit barrier at the lateral boundary. Due to 
a scaling argument that will be clear soon it is enough to consider a very special case. 

5.1. An explicit construction of a supersolution at the boundary. We dehne the func¬ 
tion 

v^{x,t) ■= \x'\^ + + {2t + 1 )_, 

where M > 1 is to be chosen depending on data. We aim to show that is a weak 
supersolution in 

Q:= {(x,f) : |x'| < 1, e [0,2], te [-1,0]}. 

Simple calculations show that 

Dv~^ = (2a:', Mx~^^‘^l2), dtv+ = - 2 x{-i<t<-i/ 2 }, 

i^\+=diag(2,...,2,-Mx-3/V4), 
and moreover, since = 0 whenever i ^ j, we have 

n n 

divA{Dv~'') = '^DiAi{Dv^) = ^ D^.Ai{Dv~^)D'^jV~^ 

i —1 

M 

= 2 ^ D^,A. {Dv+) - —D^MDv+)x-^/\ 
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The first term we estimate from above using (1.9)2 second term we can ap¬ 

ply (1.9)i, since D^^An{Dv'^) = {DA{Dv^)emen)- Furthermore, if we require M > 
23/2l6(n- l)L/z/, we obtain 


div A{Dv+) < (2(n - 1)L - (5.1) 


- 8 " \Dv+\ ■ 


Now, observe that since M > 4 we also get 

11??;+1 = \/ 4 | a ;'|2 + 

in Q. On the other hand, we have 

\Dv+\ > Mx-^/^12 > 1. 


Using these estimates we obtain 


\Dv+\ 






-1/2', go-2 


ffo > 2 
<?o < 2 


and thus 

divy4(L»?;+) < _!^M“”{9«’2}-l2.-(min{go.2} + l)/2^ 

8 

The exponent of Xn is negative, so that by choosing M = M(data) large enough (recall 
that po > 1). we finally obtain 

dtV+ - div A{Dv+) > -2 + ^2-(™{9‘>’2} + l)/2^min{go.2}-l > 

8 

It is easy to see that u+ € {Q) and thus ?;+ is a (weak) supersolution in Q. 


5.2. A reduction of the oscillation in a significant case. We set ourselves now in what 
seems to be a very particular, unitary case; it will be clear soon that, up to a simple rescaling 
procedure, this will be the significant case for the proof. 


Let fi be a bounded domain and CIt := U x (—1, 0). Suppose that 0 € dfl and 
the orthonormal system where the boundary is a graph is the standard cartesian one, with 
the direction where dCl is a graph given by e„. We hence have 

dCl n {|a;'| < 1, \xn\ < 1} = graphs, with 9 : B[{0) —>■ (—1,1) and 0(0) = 0 

and fi n {|a;'| < 1, |x„| < 1} is the epigraph of 6. Let u be a weak solution to (L8)i in 

nrnQi with Qi :=5( X (-1,1) X (-1,0) cM”+\ 

such that it = in dp^r H Qi. Moreover, we suppose '0(0) = ■i2(0) = 0. Take S G (0,1) 
to be fixed later. We assume that 

the graph of 9 over B[ is contained in the cylinder B[ x (—d, S) (5.2) 

and moreover that 

osc 'ip < S and osc ■S < 1 (5.3) 

annBix(-i,i) citCiQi 

Let us take the barrier v+ built in the previous paragraph and shift it in the Cn direction 
as follows: 

vj[x' ,XnA) '■= V'^{x' ,Xn +6,1)+ 6. 

Now v'l is defined and continuous, in particular, over the parabolic closure of Ot H Qs, 
where Qs = B[x {—6, 1) x (—1,0), and there it is still a supersolution to an equation struc¬ 
turally similar to (L8)i. The aim is to prove that 'u < on dp{ClT (T Qs) by considering 
the different pieces: 
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• on X (—1,0)] n Qs we estimate 

u — v'g < Ip — S < 0 

using (5.3)i and since u"*" > 0; 

• on [{dB[ X [—15,1]) n n] X (—1,0) we have 

by (5.3)2 together with u(0) = 0 and the fact that v~^ > 1, since |a;'| = 1; 

• on X {1} X (—1, 0) we have 

u-vj<l-M<0, (5.4) 

since u < 1 as above and on {xn = 1} we have > M > 1; 

• finally, on 17 D {B[ x (—<5,1)) x {—1} we again have u/ > 1 due to the expression 
of the time-dependent part, and therefore the conclusion again follows. 

Note that the first three pieces exhaust the lateral boundary of CIt H Qs, while the fourth 
one makes up its initial boundary. Therefore, we have u < on the parabolic boundary 
of CIt n Qs and hence, by Proposition 2.1, it < v'l in CIt H Qs- Now, if 5 < 1 /2, we have 

u+ < (5^ + M(2(5)1/2 Jjj X S)) n 17] X (-5,0). 

Therefore, if we choose <5 small enough, depending only on M and so ultimately on data, 
such that (5^ + M{2S)^^‘^ < 1/4, then we have 

1 

sup U < -. 

[(S'x(-5.<5))nn]x(-<5,0) 4 

Completely analogously we may consider the subsolution v- {x', Xn,t) = —v+{x', Xn, t) 
to obtain a corresponding bound from below. All in all, we conclude with 

osc u <^. (5.5) 

[B^nO] X (—5,0) 2 

5.3. Iteration. Let Rq < min{i?o, 1} be fixed and let (5/’(a;o, to) be a cylinder not inter¬ 
secting the initial boundary, with xq G 017, tu > 0 and r < Rq. Since we are supposing 
Bo < Bn, we have that the boundary of 17 can be written as a graph in Br'. there 
exists a unitary vector e G K" such that if we set T : M" —>■ R" for the orthogonal 
transformation that maps e„ = (0,..., 0,1) into e, we have 

T“^((?17 — xo) n (i?/ X (—r,r)) = graphs 

for some 0 G with values in (—r, r). We now start from the assumption 

osc Ue < UJ. (5.6) 

We define, for j G N, the quantities 

ujj := 2~^UJ, rj+i = min {crrj ,rj+i}, xq = r 

where <t G (0,1/2) is such that cr < -^ and (2cr)®“ < 4{\/2)~^^5 (see (5.9)), with 
6 G (0,1/2) being the constant defined in the previous paragraph, and fj is such that 

osc ifj<Sujj. (5.7) 

(xo,to) 

3 

Note that this is possible, since ip is continuous so that at ujj fixed the map p 1—>■ osCg^i ip 

vanishes as p 0. We prove by induction 

osc Ug < LUj. 

OtOQ// (xo,to) 


(5.8) 
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Now (5.8)o is simply (5.6), so we suppose that (5.8)j holds and we prove (5.8)j+i, for 
j G No. Rescale Ue as follows: 


(x, t) := i Ue (xo + ^^2;, to + UJj 


V2u 


;ixo,to) 


V2 


w,- 


This is a solution to an equation structurally similar to (1.8)i, see Paragraph 2.6, in partic¬ 
ular in [{B[ X (—1,1)) n n] X (—1,0), with boundary datum 


^{x,t) := - -0^X0 + t) - 

UJj V ^/2 I \ Tj / } J LJj 


and where the boundary of Q := [\/2T ^(17 — xo)/rj] fl {B[ x (—1,1)) is given by the 
graph of the function 0{x') = 9{rjx' jy/2)lrj over B[. We have 


- 1 1 
osc \D9\ = —= osc \D9\ < —= osc \D9\ < — 

^1 v2 V 2 ^Rg v2 


Now we choose Rq small enough so that the right hand side of the chain of inequalities in 
the above display is smaller than S, where 5 is the quantity fixed in the previous paragraph. 
This ensures that (5.2) is satisfied (since D9{0) = D9{0) = 0 = 9(0) = 9(0)). Since all 
the other assumptions in Paragraph 5.2 are satisfied (in particular by our choice of fj), we 
have estimate (5.5) at hand; therefore (5.8)j+i follows by our definition of and ujj+i- 
Indeed, scaling back we have 


osc 

^TriQrj 


1 

Ue < -UJj 


with Qr, := S 


Sr. 


/^(xo) X 



and by (2.10) and our definition of a, we infer 



(5.9) 


Finally, we note that the lengths of the time intervals also go to zero, that is, the cylinders 
are shrinking. Indeed, the first inequality in the above computation shows that the ratio of 
two consecutive time scales is bounded by (2a)^°/A, which is clearly strictly smaller than 
one. 


5.4. Some quantitative estimates. Let us set 

uj := 2||V'||l~ + 1, 

fix a radius r < Rq, and take a point (xo, to) € such that Qj^ax{i,<.j 2 / 9 o-i}r(a^o, ^o) 

does not intersect the initial boundary. Clearly (5.6) holds by the maximum principle. Now 
we recall that ijj has the modulus of continuity uj^p: 

\tp(x,t) - 'ilj(y,s)\ < w^(distpar.G((a;,f),(j/,s))) 
for all (x,t), (y,s) € dp^lx- Since 

Q-p^(xo,to) C Bp(xo) X (to - 

we have 

Q‘^^{xo,to) c Q%p(xo,to), 

with := max{l, Qj'2/go 1} > 1. Thus we see that if we want (5.7) satisfied, 

it is enough to require 

Ajfj < 
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SO that Q^.f .{xo,to) does not intersect the initial boundary, and 
< SAj^ujj 

by the concavity of oj^{-). At this point we have (5.8) at our disposal, and this will be used 
noting that in particular we have 

Qr^{xo,to) 3 Q^./B.{xo,to) 

with := max{l, > 1. Hence, for (x, f) G VIt r\ Q^^g^ixo,to) 

fixed we find the largest j G No such that 

< distpar.G((a;,f),(xo,fo)) < 

Bj+I Bj 

Note that this is possible, since clearly rj/Bj < rj —>■ 0 as j oo. At this point 

|Me(x,f) - Ue(xo,fo)| < OSC Ue<2~^UJ. 

QrriQrj (xo,to) 

Let {(rj+i/Bj+i, 2“-^ w)}jgNo a sequence of points in and call u;„ the smallest con¬ 
cave function such that w„(rj_|_i/f3j+i) > 2~^uj-, note that is a modulus of continuity. 
For instance, one can take the piecewise linear interpolation of the sequence {{xj, 
given by Xj = maxfe>j_|_i Xk/Bk, yj = 2~^uj, which is component-wise decreasing as j 
increases. This finally leads to 

\ue{x,t) - Ue{xoAo)\ < < u!u{rj+i/Bj+i) < (distpar,G((a;,f), (a^o,fo))), 

(5.10) 

and this holds for (x, f) G nTr\Qf/Bg{xo,to)- In fact, it also holds for points (x,f)outside 
Qr/Bo (xo,fo)> since then we have distpar,G((a;, f), (xq, fo)) > xIBq and thus 

\Ue{x,t) - Ue{xo,to)\ < < cdistpar,G((a;, f), {xoAo)) 

r Bq 

by the maximum principle. Note that the modulus of continuity lOu at this point depends 

on data, ||' 0 ||ioo,u;^ but also on r. 

If now Ip is 7-H61der continuous with respect to the G-parabolic metric, then we see 
that it is enough to take rj = c{data,uj,j)2~'^^r for some 77 = 17 ( 51 , 7 ). This yields 
that the numbers rj can be written as fj^r for some fj G (0,1). Now the Holder continuity 
follows, for instance, similarly to [14, Chapter III, Lemma 3.1]. 

5.5. Continuity at the initial boundary. We begin by modifying the barrier built in Para¬ 
graph 5.1 to meet the different situations at the initial boundary. We start by considering the 
case where, before rescaling, we have a solution in a cylinder Br(xo) x (0,uj^/G(uj/r)), 
with Br(xo) C fl, equal to Pp over Br{xo) x {0}; that is, the true case of initial bound¬ 
ary continuity. Later on we shall face the “corner case”, that is the case of cylinders 
Br{xo) X (OjW^/G)^/!-)) with xo G dVt. 

After rescaling, one sees that it is enough to build a supersolution in Q ;= i?ix(0,1). In 
this case the explicit expression is simply u“*'(x, f) := |x|^/^. We then have G (Q) 
and is a supersolution to (L 8 )i in Q. Moreover, if we further suppose that 

OSC 'ip < S, oscu<l, u(0,0) ='0(0,0) = 0 
Bix{0 } Q 

for some <5 G (0,1), it is easy to see that 'll < u’*' -f b on dpQ. Indeed on dBi x (0,1) we 
have = 1 but u < 1 , while on Hi x {0} we have u = 'ip < S and v~^ > 0. Therefore 
we can deduce by Proposition 2.1 that u < x"*" -f b in Q. Now the proof goes on similarly 
as in Paragraphs 5.2 to 5.4, with possibly new constants S and Rq. 

For the “corner situation”, we are lead to consider a solution in a domain of the type 
Q := B[ X (—1,1) X (0,1); the supersolution in this case is U+ = |x'p + Mxl/'^, with 
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M as in Paragraph 5.1. The fact that the function is a supersolution follows plainly from 
(5.1). Assuming now that the boundary graph 6 over B[ takes values in {—S, 5) and 

osc tp < S, oscu < 1, 12 ( 0 , 0 ) ='!/5(0, 0 ) = 0 , 

OpOrne C 

we have 22 < u/ in dpi^T H Q, since u = Pp < 5 there; on the remaining part of the 
parabolic boundary of Q we use the fact that v'l is larger than one, as in (5.4). Again, now 
the proof is similar as above. 

In both cases, a scaling and iteration procedure like the one used in Paragraph 5.4 allows 
us to prove the reduction of oscillation in a sequence of nested cylinders of the type {ft fl 
Brj {xq)) X (0, ijj'jG{u;j/rj)), with xq £ O. This leads to 

\ue{x,t) - UE(a;o,0)| < Wu(distpar.G((a:,i), (a;o,0))). (5.11) 

for every {x, t) £ fir- Moreover, opportune statements similar to above still hold in the 
case Ip is Holder continuous. 

At this point we call R the smallest value of Rq coming from the three different cases, 
ultimately a constant depending on data and dft. Choose r = R/2. Now r is a constant 
depending only on data and dft. By combining the boundary estimates (5.10) and (5.11) 
with the interior estimate in Remark 2 we obtain 


\Ue{x,t) -Ue{y,s)\ < UJu{distp^r,G{{x,t), {y, s))) (5.12) 

for every {x,t), {y, s) £ ftr^, where uJu depends on data, e, Ci,u)^, ||^||l=o, OH. Indeed, 
if one of the points is in dpflr, then (5.12) is either (5.10) or (5.11). In the case where both 
{x, t), {y, s) £ fix we consider two different cases. Either the mutual distance of {x, t) 
and {y, s) is small compared to their distance to the boundary, in which case we use the 
interior estimate, or otherwise we can again use the boundary estimates. 

Let us make this rigorous. Denotes = {x,t),w = (y, s). If distpar,G(-Zj tf) > 1, we 
are done by the maximum principle. Note now that if distpar,G(-z, w) < I, we have 

distpar,G(2, < distpar(2:, W) < distpar.G (-Z, • 


Observe that (4.15) can be written in terms of the parabolic G-distance as follows: 

\Ue{z) - Me(w)| < C distpar.G(2, 


where 

:= min{distpar.G(-2:, 9pHT),distpar.G(w,9pOT), l}. 

If now distpar,G(' 2 ^i w) < [df 27 max{i,gi/ 2 }^ mutual distance of z and w is 

small compared to their distance to the boundary, then we have 

\Ue{z) - Ue(w)| < cdistpar,G( 2 :, 

On the other hand, when distpar,G('Z; it') > [d^u)] there exists a cylinder 

Qp (xo,fo) ^ z,w with p = 2 distpar,G(' 2 ) such that either xq £ dfl or 

the bottom of Q^(xo,to) touches the initial boundary. Now using triangle inequality and 
the boundary estimates yields 

\Ue{z) - Ue(w;)| < 2a;„(p) < 4Wu(distpar,G(-2, 


Finally, we take the largest modulus of continuity for which all the conditions proved 
above are satisfied, and this proves (5.12). The proof in the Holder case is similar, since in 
this case we can quantify all the moduli. 
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6. Conclusion 


Call Uj := Ue for e = 1/j, j £ N, and similarly Aj,gj,(j)j. From the results of the 
preceding section, that is, from the equi-boundedness of the sequence following 

from the maximum principle Corollary 2.2 and the global equi-continuity coming from the 
results of Sections 4 and 5, using Ascoli-Arzela theorem we see that Uj u uniformly 
in for some u £ Now all we have left to prove is that u is a weak 

solution to (1.8) 1 , which follows easily from the next proposition. 

Proposition 6.1. Let Uj £ be the solutions to (4.3) defined above. 

Suppose there exists a function u such that Uj —> u almost everywhere in LIt- Then 
Duj —> Du almost everywhere. 

Proof. Take /C (s LIt and choose a cutoff function ip £ Cf^{LlT) such that 0 < (/? < 
I, p = 1 in /C, and \\dtp\\L’^(nT)A\Dp\\L'=°{nT) < c for some c > 1 depending on 
dist(/C, dpLlT). Let fc £ N and test the weak formulations of Uj and Uk with rj = Wj^kT^ 
where Wj^k ■= Uj — Uk. This choice can be justified by standard methods such as Steklov 
averages. By subtracting we obtain 


0 = - 


Wj^kdt{wj^kp)dxdt+ / {Aj{Duj) — Ak{Duk),D{wj^kT’)) dxdt 
' J Or 

f.dtp dx dt + / (^Aj{Duj) — Ak{Duk), Dp)wj^k dx dt 
^ Or 


+ 


(Aj{Duj) — A{Duj), Dwj^k)^ dx dt 

(^A{Duj) — A{Duk)T Dwj^k)f dx dt 

Or 

■ / {A{Duk) - Ak{Duk),Dwj^k)‘pdxdt =: I 

J Or 


+ // + III + IV+ v. 


Since ||DMjj|j;,oo(Q^) < c uniformly with respect to j by Proposition 4.1 , we also have 

ll-4j(L'wi)llL“(aT) < cgj{\\Duj\\Lo^^Q^)) < c. 

Thus, by the definition of Aj 

\I + II + III + V\< c\\uj -Mfe||L2(o^) +c||((/>j *A){Duj) - A{Duj)\\l 2 (^q.^^ 

+ c||.4(i:>Mfe) - {(j)k * A){Duk)\\L^{Q.j.) +c{l/j + 1/k). 

The first term on the right-hand side tends to zero as j, fc —oo by Lebesgue’s dominated 
convergence theorem and the second and third by the properties of mollifiers; the last one 
is obvious. On the other hand, by (2.13) 


IV > c 


[ \VADu,) 

JK 


— Vg{Duk)'^ dxdt. 


Thus 


c / \Vg{Duj)-Vg{Duk)\^dxdt<IV <\I + II + III + V\^0 

JK 

as j,k —>■ oo. We have shown that the sequence {Vg{Duj)}j^fi is Cauchy in L^{1C) 
and therefore there exists a function w £ L^iJC) such that Vg{Duj) —>■ tu in L^iJC) as 
j ^ oo. This implies that there exists a (nonrelabeled) subsequence Vg{Duj) converging 
to w almost everywhere in 1C. Now the fact that Vg has a continuous inverse yields 


Dug = Vg-^iVg{DUg)) 


Vg {w) =: V 


almost everywhere in /C. 
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Now, since Uj —?> u almost everywhere in /C, we have for any </> G (/C) that 


uD(j) dx dt = lim 


UjD<j)dxdt = — lim / Duj(j)dxdt = — / v(j)dxdt 

Jk Jk 


by Lebesgue’s dominated convergence theorem and the definition of weak gradient, show¬ 
ing that V = Du. Thus, we have Duj —^ Du almost everywhere in /C for any /C g Dt, 
which implies that Duj —^ Du almost everywhere in Dt- □ 


To conclude, (1.10) follows from (4.9) simply using the local almost everywhere con¬ 
vergence of Due, and (1.11) follows from (5.12) using the global uniform convergence of 

Ue- 


6. 1 . Weakening the assumptions. As mentioned in Remark 1, in this paragraph we show 
how to modify the proofs of the paper in order to obtain Theorems 1.1 and 1.2 for vector 
fields satisfying the weaker assumptions (1.12). 

We observe that assumptions (1.9) are only used in order to have the analogous proper¬ 
ties for the regularized vector field Ae defined in (4.1). Moreover, (2. 1 1) and (2.12) trivially 
hold by taking ^2 = 0 in (1.12). Thus, it suffices to show that under the assumptions (1.12) 
we still have (1.9) for Ae with g replaced by ge defined in (4.2). 

We shall focus only on the convolution part of the vector field Ae, since for the part 
involving the nondegenerate pi-Laplacian the corresponding estimates are classic and easy 
to verify. Therefore, we only need to prove (1.9) with A replaced hy (jje * A and g{s) 
replaced by Using (1.12) we have 

{D{4>e * A){0>^, A) = lim i [ (.A(C + hX-'g)- .A(^ - p), \)4>e{r]) dg 

h^o h Jb^{o) 

and 


\DA4>e*Am\=}^^ \ _ 


/ (.A(5 -I- hij 

> 3 ,( 0 ) 

f 9{\^-£v\) 
jBm \^-£v\ 


-g)-A{^-g))(l>e{g)dg 

4){g) dg. 


Hence, if we can show that 



g(IC-g^7l) 

\i-£v\ 


(t){g) dg 


g(iei + g) 

1^1+£ 


independently of e, we are done. 

Consider first the case |^| > 2e. This implies — eg] > |(|^| + e), and thus — 
1^1 -f e so that (6.1) holds. On the other hand, if |^| < 2£, we have 


( 6 . 1 ) 

£g\ Ri 


and 


JBiiO) — 




Bi(0) JBaii/e 


< c(n) sup 4> 

Bi( 0 ) 


g(lel+g) 

1^1+£ 



g(|g-gg|) 

\^-£v\ 


4>{g) dg > 


L 


Si/2(0)\Si/4(C/£) 


g(l^-£??l) 

|C-e?7l 


(j){g) dg 


- >l-Si/2(0) \Hi/4(^/e)| 

- 01 / 2 ( 0 ) 


g(£/4) 
1^1+£ 
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> c(n, Qi) inf 

Bi/2(0) 


gm+e) 
ICI +£ 


Note that we can assume without loss of generality that sup^^^Q^ 4> < c and inf 51 ^ 2 ( 0 ) 4> ^ 
1 /c for some c = c(n) > 0 . 
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